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Abstract 

We study the fermionic sector of a probe D8-brane in the supergravity back- 
ground made of D4-branes compactified on a circle with supersymmetry bro- 
ken explicitly by the boundary conditions. At low energies the dual field the- 
ory is effectively four-dimensional and has proved surprisingly successful in re- 
covering qualitative and quantitative properties of QCD. We investigate fluctua- 
tions of the fermionic fields on the probe D8-brane and interpret these as mesinos 
(fermionic superpartners of mesons). We demonstrate that the masses of these 
modes are comparable to meson masses and show that their interactions with 
ordinary mesons are not suppressed. 



1 Introduction 



Since the discovery of the AdS/CFT correspondence [U El G], there has been a 
lot of interest in trying to find a gravity dual of QCD. Important progress was 
made with the realisation [3] that flavour degrees of freedom can be incorporated 
by adding Nf D-branes in the probe approximation where their back-reaction is 
ignored. This is an approximation, which is valid only if the number of flavours 
is much less than the number of colours, Nf <C N c . 

Chiral symmetry breaking in holography was studied in a D3/D7-brane model 
in [5]. An interesting D4/D6-brane model was presented in ref. [6] where N c D4- 
branes compactified on a circle make a supergravity background, in which Nf 
flavour branes are put in as probes. On scales smaller than the Kaluza-Klein 
scale associated with the compact direction, the dual gauge theory is effectively 
four-dimensional and exhibits many of the qualitative features of QCD. 

Based on this model, Sakai and Sugimoto presented a similar D4/D8-brane 
model [7J [8], which comes surprisingly close to recover the meson spectrum as 
found experimentally. In particular, the spectrum contains a massless pion, as- 
sociated with a spontaneously broken chiral symmetry. Moreover, the effective 
theory for the pions was found to contain the Skyrme model. It has been argued 
that baryons can be constructed in the supergravity description by introducing 
a baryon vertex [9], which in the D4/D8 context is a D4-brane wrapped on a 
four-sphere. This object is interpreted as the Skyrmion on the probe D8-brane. 
Several investigations of the baryons from this point of view have been performed 
recently in e.g. refs. [El QH QH CGS Q31 13 [L6] . 

As a holographic description of QCD there are some fundamental difficulties 
with the D4/D8 model. One problem is that on scales small compared to the 
radius of the compact direction (i.e. in the UV) the gauge theory is inherently 
five-dimensional. Moreover, the mass scale of the scalar and vector mesons is 
the same as the mass scale of the Kaluza-Klein modes arising from the 
compact direction. That is, the model predicts a tower of non-observed Kaluza- 
Klein modes with similar masses to the observed mesons. Another problem is that 
the D4/D8-brane configuration is non-super symmetric and unstable, the unstable 
modes being the scalar fields corresponding to the radius of the circle and to the 
asymptotic brane-anti-brane separation. However, in the supergravity limit these 
modes become non-normalisable and as such cause no worry. 

The addressed problems in the D4/D8-brane holographic model of QCD might 
well be related to the fact that we are only working in the probe approximation. 
Ideally, the flavour D8-brane back-reaction should be taken into account, but how 
to do this is only partly understood at the moment. For some attempts to go 
beyond the probe approximation, see e.g. refs. [T7 1 [TH| [TTfl I2U1 I2T]. 

In this paper, we investigate an implication to the Sakai-Sugimoto model of 
string theory having a supersymmetric spectrum. Although the supersymmetry 
of the D4 background is explicitly broken by the boundary conditions, there is 
still a supersymmetric spectrum. In particular, the D8-brane probe action has 
a fermionic part in addition to the bosonic Dirac-Born-Infeld and Chern-Simons 
parts. The fluctuations of these fermionic modes should be interpreted as the 
supersymmetric partners of the mesons, so we refer to them as mesinos. Being 
aware of their presence in the full spectrum, a natural question is therefore what 
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their masses are. Since mesinos are not found experimentally, nor in the QCD, 
these masses ought to be very large. However, due to the absence of different scales 
it seems unlikely that this is so. A priori, we therefore expect mesino masses to 
be comparable to the meson masses, of order M^k- We compute these masses 
explicitly and conclude that the expectation is true: Mesinos appear on the same 
scale as mesons. 

The question then arises how this "unwanted" (from the phenomenological 
point of view) mesino sector affects the "good" meson sector. To get a partial 
answer to this question, we investigate mesino interaction terms, in particular, 
we compute trilinear meson-mesino-mesino couplings. One could be tempted to 
speculate that these interactions are somehow suppressed and that the mesino 
sector therefore would not interact with the meson sector. However, we shall find 
that this is not the case. The mesinos do indeed interact with the mesons. 

The paper is organised as follows. Section 1 contains a review of the Sakai- 
Sugimoto model, introducing the notation used throughout the paper. Moreover, 
the spectrum of the mesons in this model is re-derived. Section 2 displays the 
computation of the spectrum of the fermionic fluctuations of the D8-brane, which 
are identified as mesinos. Section 3 describes how these mesinos interact with 
bosonic fields. The effective four-dimensional interaction vertices are given ex- 
plicitly. Section 4 is dedicated to a discussion of our results and concludes the 
paper. 



1.1 The Sakai-Sugimoto Model 

The Sakai-Sugimoto model [7] is a D4/D8-brane configuration with N c D4-branes 
wrapping a circle and Nj flavour D8-D8-branes. Schematically, the brane config- 
uration is as follows: 



D4 

D8-D8 



0123456789 



The 4-direction is the circle, and anti-periodic boundary conditions for the fermions 
on this circle break the supersymmetry explicitly. 

This D4-brane supergravity background is given by [9] 



3 3 

ds 2 = (I) * (v^dxV + f(U)dr 2 ) + 5 (j*L + U 2 dnj 
r An 3R3 ,/m 1 U kk 



(i) 



where //, v = 0, . . . , 3 and r is the compact direction, which the D4-branes wrap. 
V4 = ^p- is the volume, 64 is the volume form and dSl\ is the line element of a unit 
S 4 , respectively. R and Ukk are constant parameters and R 3 = irg s N c l 3 where 



g s is the string coupling, l s = yet' is the string length and N c is the number of 
D4-branes. The coordinate U satisfies U > U^k, and to avoid a conical singularity, 
the coordinate r must have a period 

4-7T t _! 2-7T 3 1 -i 
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where we have defined the Kaluza-Klein mass Mf-k- Below this energy scale, the 
field theory is effectively four-dimensional. In the limit Ukk 0, the r direction 
uncompactifies and we recover the background built out of flat D4-branes. 

We are going to study a particular embedding of the D8-D8-branes, in which 
the D8 and D8 are antipodal in the (U, r) coordinates and smoothly join together 
at U = Ukk- For this purpose, it is useful to make the change of coordinates 
(U, r) — > (y,z), defined through 

r 2 2-7T 

U kk dT 

y = r cos 9, z = r sin 9. (4) 

In the (y, z) coordinates, the probe D8-D8-branes are extended along the z di- 
rection and located at y = 0. This configuration corresponds to massless quarks 
in the gauge field theory and was demonstrated in the original paper [7j to be a 
stable solution of the probe action. 

1.2 Meson spectroscopy 

In this section, we review the computation of the meson spectrum in the D4/D8- 
model. 

The bosonic part of the action for a D-brane probe is 

S = Sdbi + Scs + Sf, (5) 
where the Abelian Dirac-Born-Infeld action is 

Sdbi = T J d p+1 £ - det(g a p + J^p), (6) 

where g a p = d a X M dpX N Gmn is the induced metric, T a p = 2na'(F a p + B a p), 
and B a/ 3 = d a X M dpX N Bmn is the induced B-field and F Q/ g is the world-volume 
gauge field. The constant T is related to the string length l s by T = ((2-7r) 8 /^) -1 . 

The Chern-Simons part Scs is important e.g. to achieve the correct chiral 
anomaly of QCD as well as the Wess-Zumino-Witten term in the chiral La- 
grangian, but will not be relevant for us. Our focus is on the fermionic part, 
Sf, which we will return to after a very brief summary of meson spectroscopy. 

Recall that the gauge sector of the field theory is associated with open strings 
ending on the D4-branes, giving rise to a U(N C ) colour gauge symmetry, where 
iV c is the number of D4-branes. Introducing a probe brane means introducing a 
flavour sector associated with strings stretched between the D4 and the D8 or D8 
branes, giving rise to a U(Nf) global flavour symmetry, where Nf is the number 
of probe branes. 

The basic idea of the holographic description of mesons in this setup is that 
they can be identified with fluctuations of the probe D-branes. This works straight- 
forwardly for scalar and vector mesons, which are associated with fluctuations of 
the probe D-brane collective coordinates and world-volume gauge field. For higher 
spin mesons, alternative approaches have to be applied. In general, the situation 
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is complicated, but ultra-high spin mesons have been studied in refs. [23] using 
a semi-classical macroscopic spinning string description. 

Returning to the low-spin case, in order to study the fluctuations it is necessary 
to fix the coordinate invariance. This is done with a static gauge, which in our 
case is 

e = X a , for 5, X 5 (0=y(0, (?) 

where the £ a are the nine-dimensional world-volume coordinates and X M are 
ten-dimensional space-time coordinates. 

The induced metric on the D8-brane in this case is 



3 „ 3 



ds 2 m = (|) 5 V»»dx»dx» + ^fdz 2 + RMcMi (8) 



where U should now be viewed as a function of z given by 

/ 2 \ 1/3 

U(z) = U kk y- + YjT) • ( 9 ) 
We will henceforth work with dimensionless quantities 

3 3 

w ^Wk fwe(i+u,2)| ' a = {^y = {i MkkR Y • (io) 



such that the induced metric g a p reads 



ds 2 m = g a(3 dx a dxP = a 2 V^r )tlu dx'"dx u + -^-V-^dw 2 + ^-^V^dnl (11) 

M kk 4 M kk 

The D8-brane degrees of freedom that remain after this gauge fixing is the 
scalar field y (fluctuating around the value y = 0) and the world-volume gauge 
fields A a . The nine-dimensional world- volume naturally splits up in a four- 
dimensional part x^, an extra dimension w and an S 4 part. Considering only 
singlets under the 50(5) of S 4 , the gauge field degrees of freedom are A^x^, w) 
and A w (x >1 , w). 

The Abelian Dirac-Born-Infeld action can be expanded up to second order in 
the field strength as 

Sdbi = T J d 9 {e-y-det( gtt(3 +^) 



pa 



M 2 

+ (2^a'f 1 ^Vr ] ^F^F uw ) + 0(F a ), (12) 



where 



is a dimensionless constant. 
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We are interested in finding a four-dimensional effective action for the gauge 
field fluctuations with a canonical normalisation. To achieve this we expand the 
gauge fields and y in terms of complete sets {p n (w)}, {q n (w)} and {p n (w)}, 

A ll (x»,w) = Y,V { ; i) Pn (w), (14) 

n 

A w {x^w)=Y J V {n \n{™), (15) 
n 

y(x»,w) = Y,U {n H^)Pn(w)- (16) 

n 

In the following we will focus on the vector fields and assume Nf = 1. The 
expansions of the field strengths then take the form 

F^(x,w) = (V^(*) - d v V^{x)) p n (w) 

n 
n 

F, w (x,w) = Y J d^ (n) (x) qn ( W )-V^(x)p l n (w). (18) 
n 

In order to have canonically normalised kinetic terms in the effective four-dimen- 
sional action, we impose the normalisation conditions 

J dw V~ 1 (w)p n {w)p m (w) =5 nm , (19) 
TM 2 f 

/ dw V 3 (w)q n (w)q m (w) =5 nm , (20) 
and choose the p n (w) to be eigenfunctions of the equation 

- Vd w (V 3 d wPn ( W )) = \ n p n {w). (21) 



This equation together with eq. (|19j) implies 

dw V 3 (w)p' n (w)p' m (w) = X n 5 nm . (22) 



o 



Comparing eqs. (I19p and (I22|) . we notice that we can choose the eigenfunctions q n 
to be proportional to p' n : 

q n =(Ml k \ n y 1/2 p' n . (23) 

The fields (x) can be absorbed into the vector meson field by the gauge trans- 
formation 

V W(x) - V W(x) + (A„M| fc )" 1/2 d^\x), (24) 
such that the action for the vector meson field V^\x) takes the canonical form 



S = d A x 



n>l 



(25) 
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for the massive vector fields V. The field strength mode expansion is given by 

f» w {x,w) =-EnVi n) (xKH = En^(^KH- (26) 

Note that the eigenfunctions p m are chosen to be either odd or even functions 

p m (-w) = (-) m+1 p m (w), (27) 

such that the odd numbered modes are vectors whilst the even ones are axial 
vectors. 

In addition to the above massive modes, there are non-normalisable solutions 
to the eq. (f2"Tj) with eigenvalue zero: 

Po(w) = C arctan(w). (28) 

This reflects the fact that there exist a function go = CV~ 3 , which is orthogonal 
to all other q n 

J dw V 3 q q n oc J dw p' n = 0. (29) 
The normalisation constant is determined easily to be 

Taking the zero mode into account, the effective action becomes 
S = ! d 4 z k^V 0) + £ ( \f$F^ + \\ n Ml k VPv^A 

J n>l ^ 



, (31) 



where the massless scalar field is interpreted as the pion. 

The masses of the vector mesons were found in ref. [7] using a numerical shoot- 
ing technique. These masses can be matched with the observed values amazingly 
well. It is surprising because we have no really good understanding of why this 
particular D-brane setup should be a good approximation of the holographic dual 
of QCD. 

So far, we have tacitly assumed that the gauge potential vanishes as w — > 
ztoo. The field strength should vanish as w — > ±oo for the effective action to be 
normalisable. Then it is always possible to choose a gauge where the potential 
vanishes for large \w\. In section l3~T1 we shall study couplings to an external gauge 
field, and for that it is useful to employ the gauge where A w = 0. At first look 
the pion field ip^ seems to be gauged away in this gauge but this is impossible. 
There exists no gauge, which sets A w = and lets A^ vanish for large at the 
same time. 

The A w = constraint can be accomplished by the gauge transformation 

A a — ► A a — d a A (32) 
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with the gauge function A given by 

oo 

A(x,w) = ^(x)Po(w) + Y,<P (n) (x) (M 2 kk \ n y 1/2 p n (w). (33) 

71=1 

After the gauge transformation eq. ([32]) . the gauge potential reads 

A,(x,w) = -d^\x) Po (w) + ^2(v^\x)- 1 — d^ n \x))p n (w), 

A w (x,w) = 0. (34) 

The contribution from the massive scalar modes can again be absorbed into 
the massive vectors by a field redefinition. The asymptotic values of the 
gauge potential for large values of |w| are 

A„(x,-oo) = C^d^(x), (35) 

A fl (x,oo) = -C^d^°\x). (36) 

Taking an external gauge field into account, the gauge potential in A w = 
gauge may be written as 

A^(X,W) = %+A^ P0 (w)+Y, v i l) P2n-l(w) + Y, a PP2n(w), (37) 

n>\ n>l 

A w {x,w) = 0, (38) 

where the fields Dp and ajf* are the massive vector and axial vector mesons, 
respectively. The zero modes are given by 

1 i 
A* = n ( A Ln ~ A Rfi ) + -rd^ir, (39) 

" In 

V M = \{A Lp + A Rii ), (40) 

where A^^, Ar^ are the asymptotic values of the external gauge field for w — > 
±00 and the pion field tp^ = in with its decay constant f w = 2/(tt C). In 
particular, the electromagnetic photon field can be extracted by setting Al^ = 
Ar^ = eQA^ M for the electric charge Q. The field strengths read 

F^x, w) = dn%(x) - duV^x) + (d^Auix) - duA^x)) p (w) 

+ E " 9uv { ;\x)) p 2n ^(i 

n>l 



+ E - d„a^(x)j p 2n (w) 

n>l 

= E F ff(^H' ( 41 ) 

n 

F^ W {X,W) = A„(x)p' {w) + E U i n) ( X )^2n-lH + E a i n) ( X )^2nH 

n>l n>l 
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where we indicated the split into vector and axial components. Equipped with 
these expansions we can aim for effective couplings of mesinos to gauge fields. 



2 Mesino spectroscopy 



As promised in the introduction, we shall now turn to the fermionic sector of 
the probe D-branes and compute the masses of these fluctuations. A similar 
calculation has been partly done in the appendix of ref. [7J. The same ideas have 
been applied also in ref. [M] for computing the masses of the meson superpartners 
in the D3/D7 setup. In that model, M = 2 supersymmetry is preserved and the 
masses can be calculated analytically. As expected, there is a match between 
fermion and boson masses. 

The fermionic part of the D8-brane action is given by [25|. [26] [27] 

S& =iT J d 9 ^-(g + ^)e-^^(l - r D8 )(r a D Q - A + L D8 )*, (43) 

where 



p _ \^~9 r (0) r ll^ ( 1 r a 1 ...o l 2 a T <r 



«2q-l"2q ' 



(44) 



Lds = ^_ ( ^ } r S E { q ^j^-i r" 1 """ 1 ^^ ■ ■ ■ %, (45) 

rg 8 = ^7^r a ,.. Q9 = -r 11 ^, (46) 

A = \T M d M 4> - l^F NPQR T Np Q R , (47) 
D M = V M - ~e*F NPQR r NP Q R r M , (48) 
V M = d M + (49) 

The Si's are the spin connections, V m = + \^ AB ^abm is the usual covari- 
ant derivative and F NPQR = rl N F p T Q T R l Latin indices M = 0,... ,9 refer to 
space-time and Greek indices a = 0, . . . , 4, 6, . . . 9 refer to the D8-brane world- 
volume. Quantities are pulled back onto the D8-brane in the usual way, e.g. 
V Q = <9 Q X M V m- Transformations between curved indices a, M and flat (lo- 
cal Lorentz frame) indices a,M_ are done with the vielbeins, so we have e.g. 

T a = g a @ dpX M E M ~T m_- 

The field ^ is a 32-component Majorana spinor in ten dimensions and 

P- = l(i-r D8 ) (50) 

is a projection operator necessary for kappa invariance. The T^s matrix is in the 
case of vanishing gauge field {T a p = 0) 

r D8 = I J^i-^rv.^r 11 = rv (51) 
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We are dealing with Majorana spinors in ten dimensions. These spinors obey 
a reality condition. For the product of two Majorana spinors with an arbitrary 
number of T matrices, we have 

-pMi...Af n ^ = / \n(r*+l)/2^pAfi...Af ^ 

so that these expression vanish if the spinors are the same and n = 1,2, 5, 6, 9, 10. 
Moreover, choosing as the n symmetry fixing condition the chirality condition 
r 11 ^ = \&, then the only non- vanishing products are 

since the chirality matrix T 11 transforms products of four and eight V matrices 
into six and two, which vanish. 

As mentioned before, setting y = is equivalent to setting the quark mass to 
zero. This is an assumption we are making throughout this paper. Since we are 
now not interested in the fluctuations of the bosonic fields, we set J- a p = in the 
calculation of the mesino spectrum and also set dy = 0. We will include the gauge 
field contributions when considering mesino interactions with vector mesons in the 
next section. We find that the contributions from the flux F A cancel in the action 
(I43p . The measure and dilaton factors together give 

T^e~* = f{V A {2^a') 2 y l aV 2 (54) 

With the definitions given above and some rather straightforward manipulations, 
we find 

S = if(V 4 (2ira') 2 )- 1 J d^dwd^V 2 £, (55) 



where 



C = HP- 



1_™ q4 „ 3, 

(56) 



5 „ 13 7 
+ M kk V*T±d w + — M kk wV—*T ± 



where m = 6, 7, 8, 9 and T— Vj^ is the Dirac operator on a unit four-sphere. 
We rescale the spinor 

$ # = v~ n/8 ^. (57) 



This rescaling removes the last term in Eq. (l56p . Moreover, it leads to an action 
where the term with the w derivative has no weighting factor V 



S = - %T — n I d^xdwdQA VP. 



|M fcfc y-sr^v^ 4 + v- 2 rnd^ + M kk r±d w 



(58) 



V A (2Tra') 2 

2.1 Spinor decomposition 

We want to decompose our spinor ^ into an S 4 part a 3 + 1-dimensional part 
tjj and a remaining two-dimensional part u, i.e. 

^ = u®ip(x,z)®x(S 4: ). (59) 
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A convenient choice for the decomposition of the 10-dimensional gamma matrices 
is 

^ = ^0^81, ji = 0,l,2,3; (60a) 
r^ = o- (B ®7®l J (60b) 
= tr,, ® 1 (8) 7, (60c) 
= a y <8 1 ® 7-, m = 6, 7, 8, 9; (60d) 

where 7 = 27-7^7-7- is the chiral matrix in 3 + 1 dimensions and 7 = 7^717-7- 
is the chiral matrix on the tangent space of 5 4 . In this decomposition, the ten- 
dimensional chirality matrix takes a particularly simple form: T 11 = a z (8 1 (8 1. 
If we further decompose the u into a z eigenstates u±, we get 

a z u± = u±, a x u±=u T , a y u± = ±iu T . (61) 

The kappa symmetry fixing condition r 11 *!' = ^ is then equivalent to u = u + = 

IT)- 

The general solution for the Dirac equation on the four sphere is well-known 

EE], 

X ±ls = i\f X ±ls ; \f = ±(2 + l); 1 = 0,1,-.., (62) 

where the quantum number s = 0, 1, . . . , d\ and d\ = 4(3 + Z)!/(3! /!) is the degen- 
eracy of the eigenvalue A, . With the above splitting and ^> = ( j \a x ®vp ® x\ the 
Lagrangian becomes 

>C~ <t I T> X f X (63) 
Normalising the fields according to 

(TIT)/ dSU (x< s )V V = (2ira') 2 V A 5 l > l '5 s ' s ', (64) 

we arrive at the action 



S = if I d 4 xdw ip 



~M kk \V-l + V- 2 ^ + M kkl d u 



tp. (65) 



2.2 Four-dimensional action in canonical form 

In order to read off the mass of the four-dimensional fluctuations, we need to 
rewrite the action (|65|) in canonical form. We do this by working with two-spinors 

*-(£)■ < 66 > 

and choose the Weyl basis where 

where = (1, — a 1 ) and <7^ = (1, cr J ). Then we expand the fermion fields in terms 
of complete sets {f+(w)} and {f-(w)}, 

M^,™) =^2$\x»)fi(w), if,-(x»,w) =52il><?\x»)f2(w), (68) 
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where the functions /? and are real eigenfunctions of the coupled first-order 
differential equations 



2A 



V-2fl( w )+d w fl{w) =M n V- l r_{w) 



3 

2A_. > 



where M n = jj^- With the normalisations 



(70) 



fjdw V- 2 (w)mw)f¥(w)=5 nm , 
the action takes the canonical form 

S = -J d 4 ^{^i n)t m^i n) + ^Ha^d^f + M n [^ n) Vi n) + ^fyi* ] } 



which is the action, written in the Weyl basis, for a set of Dirac spinors 

**-{$) 



(71) 



(72) 



with masses M n . With these Dirac spinors, the action takes the familiar form 



S = i J d 4 xj^ [| (n Vd^ (n) + M n ^^ 



(73) 



In order to compute the masses M n we have to solve the eigenvalue problem (169)) . 
We will do this numerically in the next section. 

The coupled first order differential equations (|69|) can be rewritten as two 
decoupled second order differential equations in Sturm-Liouville form, 



d w [V 2 d w fl{w)} +q±Wf±W = M 2 n V- 2 fl{wl 



where 



H ■ 1 9 9 36 



IA\/1 + -^±1 



(74) 



>0. (75) 



We therefore know that the eigenvalues M n are real and well-ordered and that 
the eigenfunctions are real and unique. At this point we note that A — > —A is 
equivalent to swapping f±—*f£, so we can therefore assume A > without loss 
of generality. 

Inspecting the equations, we note that fH(—iv) satisfies the same equation as 
f™(w). The uniqueness of the eigenfunctions then implies that they are the same 
up to a sign, i.e. we have 

f^- w ) = c^fl{w), C W = ±1. (76) 

For the special case M = 0, the equations (fB1?|) can be solved analytically, 
giving 

f ± = C±(w + yV + l^f . (77) 
However, this is a non-normalisable solution. 
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2.3 Numerical results 



The asymptotic behaviour of the eigenfunctions is found by investigating the dif- 
ferential equations (|69|) for uu — » oo and yields 

(2A 1 4A — 1 2A 
A W 3 3 — B ~ — W 3 
3M , 

, ( a 4A + 1 _2A 2A_l\ V ^ 

/_ ~ -A = 10 3 + 5 W 3 3 

V 3M / 

where A and .B are undetermined constants. As we have already noted, it is 
sufficient to consider the A = (I + 2) > 2 case. Then we have 

f + f + ~&w*, /_/_ ~ B^i^" 1 ) (79) 

This expression is not normalisable according to equation (|7Up unless we require 
B = 0. It is this regularity condition, which produces a discrete mass spectrum. 

Let us now turn to the initial conditions. Evaluating the equations (|69p at 
w = gives 

±/r(0)-y/£(0) = M„/^(0). (80) 
From equation ([75]) we get 

/_(0) = c/+(0), /l(0) = -c/;(0). (81) 
Equations ([80l) and ([HT]) together give 

(82) 



/+(0) = /o; A(0) = (f -cM)/ ; 
/-(0) = c/ ; /i(0) = -c(f -cM)/ , 



so that the only choice we have in initial conditions is in /q and the sign of c. 

The aim when solving the eigenvalue problem numerically is to identify the 
mass eigenvalues M n corresponding to regular eigenfunctions. In doing this we 
apply a shooting method. The actual value of /o is not important and can always 
be set to /o = 1. The recipe is then to choose either c = +1 or c = - 1 and guess a 
value of M. Then we integrate the equations numerically up to very large w and 
compute the coefficient B in the asymptotic behaviour ([78]) according to 

B = lim w^ + ^f_( w ). (83) 

w— >oo 

Different M give different values for B, and as we have argued above, only those 
giving B = are regular mass eigenstates, see figured) In order to determine all 
the eigenvalues, we scan through a range of M for c = ±1 and for various / > 0. 
It turns out that the lightest masses are found for c = — 1. In general, odd (even) 
n corresponds to c = — 1 (c = +1): 

f n { _ w) = (-)»/£(„,). (84) 

Numerical values for M are presented in Table [TJ 

We can determine the parity of the eigenstates in the same way as in ref. [7J. 
Consider a transformation L : (t,x,uu) — * (t,—x,—w), which is a proper Lorentz 
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1 . 5 

(M 2 + l)B 

1 



0.5 



-0.5 



AI 



Figure 1: The coefficient B plotted as a function of M for A = 2 and c = ±1. The 
intersection points B = give the mass eigenvalues. 





n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


/ 


= 


1.60- 


2.11+ 


2.59- 


3.05+ 


3.51 


3.96+ 


/ 


= 1 


2.27~ 


2.80+ 


3.29- 


3.77+ 


4.24" 


4.70+ 


/ 


= 2 


2.94- 


3.48+ 


3.98- 


4.47+ 


4.95- 


5.41+ 


/ 


= 3 


3.61- 


4.15+ 


4.67- 


5.17+ 


5.65- 


6.13+ 



Table 1: Mesino masses in units of M^k with parity indicated. 



transformation in five dimensions. The four-dimensional interpretation is that of 
a parity transformation. Acting on a spinor, the transformation is 

ip(t, x, w) -> (-i^tptt, -x, -w), (85) 

which, using the expansion (|68p gives 



(86) 



where we have used the relation (|84p . In terms of the mesino field defined in 
equation (f72|) . we see that this has the effect 

* w ~ (?) - '->" ($) ^ <-»"<-- a » (?) ■ < 87 > 

From this result we see that even-numbered modes correspond to mesinos with 
even parity, while odd-numbered ones have odd mesino parity. The parity of the 
lowest mass eigenstates is indicated with a sign in Table [TJ 

For comparison, the massless scalar meson (the pion) has negative parity, as 
does the lightest vector meson (p), whilst the lightest scalar meson, ciq(1450), has 
positive parity. See ref. [7]. 

It is interesting to compare the mass spectrum to that of the ordinary mesons. 
The mass square ratio between the lightest massive mesino and the lightest massive 
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meson (the p meson, see ref. [7]) is 



M 2 _ l-60 2 M| fc _ 
m 2 0.67M 2 , 

In comparison, the mass square ratio of the second lightest and the lightest meson 
is 2.4 (from this D4/D8-brane model) / 2.51 (experimentally). What we have 
found is that the mesino masses are comparable to the meson masses. 



3 Mesino interactions 

We have found that the model predicts fermionic partners to the mesons, which 
we have called mesinos. Their masses are of the same order as meson masses, yet 
they have no counterpart in QCD or in experiments. This by itself presents a 
problem for this model as a holographic dual of QCD. However, there is still the 
possibility that these mesinos do not interact with other particles. If this is the 
case, we would be able to view this fermionic sector as irrelevant since it does not 
affect the meson sector. This section is devoted to the question of whether this 
is in fact the case. Specifically, we will investigate the coupling of the mesinos to 
vector mesons as well as to an external gauge field. 

As in previous sections, we shall consider only SO (5) singlet states. Thus we 
set A m = 0, m = 6, 7, 8, 9 and choose the other five components of the gauge field 
to be independent of the coordinates of the S 4 part. 

If a gauge field is switched on the physics on the D8-brane is described by 
the Dirac-Born-Infeld action, eq. ©. We are interested in trilinear Yukawa- 
like couplings containing just one boson. The commutator in a non-Abelian field 
strength necessarily introduces a second boson. This would only be relevant if we 
would also consider two boson - two fermion couplings. For this reason there is 
no loss of generality by considering an Abelian gauge field. 



3.1 Trilinear Couplings 

Trilinear Yukawa-like couplings arise from terms of first order in the field strength 
in the fermionic D8-brane action (|43|) , Both Tjjg and Ld$ contain a piece linear 
in T 

F m = THl-h ll T a ^ aP + 0(F 2 )), (89) 

L m = -T^-T^FjDp + OiF 2 ) (90) 

Inserting the expressions into the action, we obtain the Yukawa-like coupling terms 
as 

-4 / d 9 £V=5e-^(i - v^T n v a ^D^. (91) 

This vertex is simplified by imposing the Majorana and Weyl condition, i.e. just 
keep terms with derivatives or three or seven T matrices. The first piece just 
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contributes 



S. 



(i) 



-i? / d^^e-^r^^na^ 



(92) 



Decomposing the spinor in the same way as in the determination of the spectrum 
eq. (|59p . we see that this term vanishes. Expanding the second piece gives the 
following 



S. 



(2) 



-iT I d 9 £ V 
1 



Vfr I (1 _ r ^ 



{T^rfPdp + M 2 kk V 4 T»T^ w 5 ww d w + M kk V 2 T^T wllV ^d p ) * 



8a 3 



16a 2 



^6789 ( y-^iw^ + 2M kk y l / 2 v^LT aw \ r 



^56789 



(93) 



From now on, the raising and lowering of indices is understood to be performed by 
rjuu and S ww . We also rescale \P — ► V~ 13 ^ 8 ^/. This has no effect except changing 
the overall weighting function. The term arising from the w derivative can be 
dropped since it comes with either one or two T and does not contribute due to 
the Majorana property of the spinors: 



5, 



a 2 y 4 (2vra / ) 



d xdwdfl^ 



frl (1 - Pi) 



W/l'- 



We would also like to add a total derivative term to the interaction 



(94) 



5, 



TM, 



kk 



2 a 2 V 4 (2Tra') 2 



™lk 



d 4 xdwdf} 4 d w [■$- (i - r§) yir/i^^ 



+f 1 (i - r^) ylr^^^ 



(95) 



in order to make the action more symmetric regarding the w derivative. The 
terms where the derivative hits the function V or F^ w vanish due to the Majorana 
condition. 
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Now we employ the splitting (|59p we used for the spectrum, perform the inte- 
gration over the four-sphere and use the normalisation for the u and x spinors 



5, 



" Z 2^2 ' d xdw 



4 



-w 



4 
M kk 



We also re-express the flux 



3i? 3 /3 a 



6789 



i\-4 2 M kk , 



(96) 



(97) 



g s \2M kk J g s a" 

As next step we go to the Weyl basis according to eqs. ([66]) and ([67]) and 
expand the fermions and bosons in terms of complete sets. The effective four- 
dimensional action is then 



m,n,p 



+t m ,n, P Ml k Fff ^ha^f + {-f^ia^) 
+lm,n, P M kk Fff (^ft^M + {-f^d^f) (98) 
+s m , n , p M kk F^ + (-)Vi B) W?' 



where & = m + n + p-|-l determines the parity properties of the effective vertex 
and Mj~^ = 27ra / (aT) _1 / 2 /2 is a dimensionful quantity setting the scale for the 
interactions. The purely numerical coefficients are given by 



Jm,n,p 

l"m,n,p 
Sm,n,p 



dw V- 7 / 2 p m (w)f^w)fl( Wj , 



dw \wy m (nu p - tin + ^- 3/2 m 



(99) 
(100) 

dw V- 3 / 2 p' m (w)f^(w)f p _(w), (101) 

dw i (w V - g / 2 + iV- 3 ) Pm (w)nHf p ~H- (102) 



The coefficients have been simplified by using the parity properties of the eigen- 
functions f± eq. ([76]) and p m eq. ([27]) . 

Written in terms of the Dirac spinors (|72p , this gives 



+M kk F$ (l m , n J {n \-7) k d^ +t m!n , p M k J^Y(l) k ^ ip) 



(103) 
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77? 


,n,p 


Jm,n,p 


L m,n,p 


7 

l m,n,p 


°m,n,p 


1 


i 


i 
i 


0.42 


0-1-0 026? 





0- 11? 


1 

1 5 


1 

3 


2 


073 


n noq — 064? 


1 1 


_n 03Q — 042? 


1. 


2, 


2 


0.26 


0.0-0.006H 


0.0 


0.0-0.051? 


2. 


1, 


1 


-0.068 


0.0 + 0.17? 


-0.33 


-0.028 


2. 


1, 


2 


0.16 


0.24 + 0.052? 


0.13 


0.011 - 0.049? 


2. 


2. 


2 


0.029 


0.0 + 0.056?; 


-0.040 


-0.041 



Table 2: The numerical values of the coupling coefficients for the interactions of mesinos 
with (axial) vector mesons. 



Since this action is effective, the non-renormalisability of the first three terms is 
not a problem. The last term is a marginal term, which minimally couples the 
(axial) vector mesons to the (axial) vector current depending on the value of k. 

Due to the properties eq. (p7|) and eq. ([84]) the coupling constants satisfy the 
following symmetry properties 



jmpn — jmnpi Impn — ( ) I 



mnpi 



tmpn ~ tmnpi s mpn ~ ( ) s mnp- 

The integrals in these expressions can be computed numerically and generi- 
cally produce numbers of order 1. Table [2] lists these numerical values for some 
interactions. What is important for our purpose is not the exact values of these 
coefficients, but the fact that they are non-zero and not very small. 

The value of the coupling constants are determined by these coefficients as well 
as by the factors Mj„£ and M^. The string coupling g s and the constant a can 
be expressed in terms of gauge theory quantities gym and A^ c [7]: 



m 



the regime 1 <C g Y M N c < 9 Y m M giving 



Mi nt = ^{ 9YM N c fN c ML- (105) 

Notice that the only adjustable scale here is M^k- The coupling constants are 
therefore given by a fixed numerical factor times the appropriate power of M^k- 
Note that according to eq. (|105f) the interaction scale Mj„f is proportional to the 

— 1/2 

fourth root of the number of colours so that the interaction terms scale with N c 

To fix the numerical factors, we can use the condition that the rho meson mass 
and the pion decay constant should match experimental values. The experimental 
rho meson mass is 776 MeV, which gives M^k = 949 MeV. The pion decay 
constant f n is related to the normalisation constant C in eq. (|30f) . giving 

ft = ^(9 YM N c )N c M 2 kk . (106) 
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Using N c = 3 and matching this to the experimental value f n = 92 MeV gives 
A = gy M N c = 8.2. With these values we find M int = 173 MeV = 0.18 M kk . The 
't Hooft coupling satisfies the requirement to be larger than one, which is necessary 
for using classical supergravity only. This restricts to small curvatures as usual in 
the AdS/CFT correspondence. On the other hand, the requirement A <C gy-M ls 
not met. It arises from the requirement that the string coupling g s = e^ ought to 
be small, which can only be satisfied up to a critical radius U cr it [6]. The main 
reason for this failure is that we are using N c = 3 while the N c — > oo and gyu —* 
in the 't Hooft limit. If we go to the 't Hooft limit with finite A the interaction 
scale goes off to infinity and the interactions are suppressed. 

Since the cross section for a particular interaction is proportional to the cou- 
pling constant and we have found the coupling constants to be non-zero and 
non-suppressed, we conclude that these interaction are not suppressed. 

Other interesting quantities are the couplings of the mesinos to the external 
electromagnetic field and to the pion. We use the zero modes in the expansions 
eqs.([HJ and (pEZj) . The mesinos couple to the external electromagnetic gauge field 
through 



iMr^ / d 4 x (d^x) - d u V u (x)) (jv,n, P ^ n h^T +P d lJ ^ ) 



n.p 

with the coefficients given by 

JV,n,p = JdwV-WfWfM, 



SV,n,p 



\f dw (wV-W + iV-^f^f® 



(107) 

(108) 
(109) 



The mesino interaction with the pion field is described by 



2 

int 



n,p 



d M 7r {M kk t^ n ^ir(l) n+p+1 ^ (p) 
+Wl (n) (-7) n+p+1 c>^ (p) 



with the coefficients given by 

L,n,p = jdwV-^f^f^\ 



(110) 



111) 



u Tr,n,p 



\ J dw (/f/; (p) - /; (n) /i p) + ^- 3 / 2 /i n) /i p) ) v-v\ (H2) 



The numerical values for the lowest mesino excitations are displayed in Table El 
As in the case of the massive modes the couplings are less than unity but still 
significantly large such that they are not suppressed and have to be taken into 
account. Mesinos couple in the Sakai-Sugimoto model to the photon and the 

— 1/2 

pion. The couplings fall off as N c and 1/N C for the photon and for the pion, 
respectively. 
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n,p 


J\P Ti.p 


<^V 71, 7) 






1,1 


0.79 


-0.20« 


0.2H 


-0.38 


1,2 


0.096 


0.089 - 0.080i 


0.031* 


0.15 


2,2 


0.58 


0.13i 


0.14i 


-0.17 



Table 3: The numerical values of the coupling coefficients for the coupling of the mesinos 
to the electromagnetic field and the pion, respectively. 

4 Discussion 

In this paper, we have investigated the fermionic sector of the flavour brane in the 
Sakai-Sugimoto model in the case of massless quarks and only one quark flavour, 
Mf = 1. We found that there is a spectrum of mesinos (fermionic mesons) with 
masses on the same scale as the bosonic mesons, set by the compactification scale 
Mkk- The presence of this fermionic sector itself is a problem for this model as a 
holographic descriptions of QCD as it is not part of QCD nor seen experimentally. 
However, one could have hoped that the fermionic sector drops out in the low 
energy physics because of very high masses. Our result shows that this is not the 
case - mesinos appear on the same energy scale as mesons. This follows from the 
fact that the supersymmetry breaking scale is the same as the compactification 
scale or meson mass scale, M^k- 

Having found mesinos of similar mass as the mesons, we asked whether their 
interaction with the mesons is suppressed. If that had been the case, we could 
have argued that the fermionic mesinos represent harmless junk in the model, 
not affecting the meson physics. We found explicitly that this is not so for finite 
N c . The interactions between the mesinos and the mesons are not suppressed and 
cannot be ignored. The existence of mesinos that interact with the mesons is a 
serious problem for the Sakai-Sugimoto model. This should in fact not come as a 
surprise, as there is only one energy scale, Mkk, in the model. 

Most treatments of the D4/D8-brane model consider only the bosonic part of 
the probe brane action and therefore do not encounter the mesinos. If the aim with 
the holographic model simply were to reproduce QCD-like physics that would be 
justified. However, one of the main motivations for studying this model is that it is 
a superstring theory model where the D-branes are equipped with a supermultiplet 
of fields. We have in principle no choice but to include the fermionic sector of the 
D-brane probe action. 

The explicit breaking of the supersymmetry due to boundary conditions on 
the D4-branes does not affect this point. The number of degrees of freedom in the 
particle spectrum is still that of a super symmetric theory but the mass degeneracy 
is lifted as we saw in our calculation. 

As said already, the problems exhibited in this paper stem from the fact that 
the supersymmetry breaking scale is the same as the meson mass scale. One solu- 
tion to it could therefore be to disassociate these scales, i.e. have a supersymmetry 
breaking scale different from Mkk • Then we could achieve a separation of mesino 
masses and meson masses. 
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